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ABSTRACT 

We present an analysis of a Chandra observation of the massive, nearby galaxy cluster Abell 2319. A 
sharp surface brightness discontinuity — suggested by previous, lower angular resolution X-ray imaging — 
is clearly visible in the ACIS image. This ^--^300 kpc feature suggests that a major merger is taking 
place with a significant velocity component perpendicular to the line of sight. The cluster emission- 
weighted mean temperature is 11.8±0.6 keV, somewhat higher than previous temperature measurements. 
The Chandra temperature map of A2319 reveals substructure resembling that anticipated based on 
hydrodynamic simulations of cluster mergers, and shows an associated cool core not previously known. 
The map shows a separation between the intracluster medium (ICM) and galaxies of one subcluster, 
indicating a transient state in which the ICM has been stripped from the subcluster galaxies (and 
presumably the dark matter). Detailed analysis of the merger feature shows a pressure change across 
the surface brightness discontinuity by a factor of < 2.5. The higher density side of the front has a lower 
temperature, suggesting the presence of a cold front similar to those in many other merging clusters. 
The velocity of the front is roughly sonic. 

We compare bulk properties of the ICM and galaxies in A2319 to the same properties in a large 
sample of clusters as a way of gauging the effects of the major merger. Interestingly, by comparing 
A2319 to a sample of 44 clusters studied with the ROSAT PSPC we find that the X-ray luminosity, 
isophotal size, and ICM mass are consistent with the expected values for a cluster of its temperature; in 
addition, the i^-band galaxy light is consistent with the light-temperature scaling relation derived from 
a sample of ^100 clusters studied with 2MASS. Together, these results indicate either that the merger 
in A2319 has not been effective at altering the bulk properties of the cluster, or that there are large but 
correlated displacements in luminosity, isophotal size, ICM mass, galaxy light, and emission-weighted 
mean temperature in this cluster. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Abell 2319) — X-rays: 
galaxies: clusters 



1. INTRODUCTION 

Galaxy cluster mergers are highly energetic events, driv- 
ing shocks into the intracluster medium (ICM) of the col- 
, liding clusters. Flattened and asymmetric X-ray mor- 
■ phologies are signatures of recent merging (Mohr et al. 

1993), and these signatures have been used to study the 
, prevalence of merging in large samples of present-epoch 
clusters (Mohr et al. 1995; Buote & Tsai 1996). A study 
of X-ray images of a flux-limited sample of 65 clusters in- 
dicates that more than half of nearby clusters show ev- 
idence of merging (Mohr et al. 1995). Hydrodynamical 
simulations indicate that complex temperature structures 
should also be produced in these mergers; however, un- 
til relatively recently the required spectral and angular 
resolution to map this structure has not been available. 
Chandra and XMM-Newton, with their high angular res- 
olution, are well-suited for detailed studies of merger fea- 
tures in galaxy clusters (e.g., Markevitch et al. 2000; 
Vikhlinin et al. 2001; Markevitch & Vikhhnin 2001; Sun 
et al. 2002; Markevitch et al. 2002; Kempner et al. 2002; 
Maughan et al. 2003). These studies have already revealed 
that merger features observed in clusters may not indicate 
shock fronts, but rather "cold fronts," wherein the cool, 
dense cores of clusters survive through the initial impact 

^ Department of Physics, University of Illinois, 1110 West Green St, 
^ Department of Astronomy, University of Illinois, 1002 West Green 
^ Current Address: Instituto de Astrofi'sica de Andalucia, CSIC, Spj 



of the merger (Markevitch et al. 2000). In fact, it now 
appears that many well-known merger candidates contain 
these cold fronts, e.g., A2142 (Markevitch et al. 2000), 
A3667 (Vikhlinin et al. 2001), A2256 (Sun et al. 2002), 
and A85 (Kempner et al. 2002). 

Abell 2319 is a massive nearby cluster (z = 0.0564; Abell 
1958; Struble & Rood 1987). We chose to study it with 
the high resolution of Chandra because its X-ray morphol- 
ogy observed at lower resolution with the ROSAT PSPC 
shows a strong asymmetry or "centroid variation," which 
is a classic indicator of a recent merger. Our goal in this 
study is not only to better understand the merger state of 
A2319, but also to determine how the ongoing merger in 
A2319 is affecting its bulk ICM and galaxy properties. Of 
particular interest is understanding how merging — which 
has long been known to be prevalent in the cluster pop- 
ulation (Geller & Beers 1982; Dressier & Shectman 1988; 
Mohr et al. 1995) — is likely to impact attempts to use clus- 
ter surveys to study cosmology (e.g., Haiman et al. 2001; 
Randall et al. 2002; Majumdar & Mohr 2003a,b; Hu 2003). 

On the basis of galaxy spectra, Faber & Dressier (1977) 
suggested that A2319 is actually composed of two clus- 
ters superimposed along the line of sight, with the smaller 
subcluster located ~ 10' to the northwest of the main clus- 
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ter and X-ray surface brightness peak. Additional redshift 
measurements led to an estimated mean veloeity for the 
main subcluster of ~ 100 members (hereafter A2319A) 
of 15727 km s~^, and for the smaller subcluster of ~ 25 
members (hereafter A2319B) of 18636 km s~^ (Oegerle 
et al. 1995). This analysis suggests that there is a ~ 50% 
chance that the two subclusters are in fact not gravita- 
tionally bound and will not merge. 

A2319 has been extensively studied with previous X-ray 
instruments, and the inferences about the cluster dynami- 
cal state have been quite varied. Emission-weighted mean 
temperature estimates are generally in the 9-10 keV range 
(e.g., David et al. 1993; Markevitch et al. 1998; Molendi 
et al. 1999; Irwin & Bregman 2000). Markevitch (1996) 
produced a temperature map of A2319 using ASCA. These 
observations provided no evidence for a cold core region 
near the surface brightness discontinuity, though a region 
to the northwest of the brightness peak appeared to have 
a temperature ~ 1.5 keV lower than the mean. This same 
subcluster region was identified by Molendi et al. (1999) 
using BeppoSAX; it is proposed that this cool region is 
associated with subcluster A2319B. Using the ASCA tem- 
perature map, Markevitch (1996) argued that there is no 
evidence of a large-scale merger in A2319. Mohr et al. 
(1995), however, found a value for the centroid variation 
of A2319 in the Einstein IPC image that indicates an on- 
going merger. Interestingly, a combined X-ray and radio 
study of the cluster suggests that the two subclusters are 
in a premerger state (Feretti et al. 1997). This study also 
takes note of X-ray evidence for another merger in a late 
stage taking place along the northeast-southwest direction. 




Fig. 1.— ROSAT PSPC image of A2319 with Chandra obser- 
vation footprint overlaid. North is up and east is to the left in all 
images. The ACIS-I footprint is roughly 17' on a side. 

In this paper, we present a detailed X-ray study of 
A2319 based on imaging spectroscopy from Chandra 
ACIS-I, providing clear evidence for an ongoing merger 
of two major subclusters. In §2 we present the observa- 
tions. After a description of the data reduction process in 
§2.1, we present an analysis of the overall cluster spectrum 
(§2.2) and a temperature map of the cluster (§2.3). In §3 
we analyze the merger feature in detail, including quan- 



titative estimates of changes in the physical state of the 
ICM across the feature, and propose a simple dynamical 
model. This is followed in §4 by a study of how this merger 
has affected the bulk X-ray properties of the cluster; we 
examine how A2319 — a cluster in the middle of a major 
merger — behaves relative to an X-ray flux-limited sample 
of clusters in its luminosity, isophotal size, and ICM mass. 
Finally, we summarize our conclusions in §5. 

Throughout the paper we assume a ACDM cosmology 
with f^M = 0.3 and Vtjy = 0.7, and take the Hubble pa- 
rameter to be Ho = TOhro km s^^ Mpc~^. 

2. OBSERVATION 

A2319 was observed with Chandra on 15 March 2002 for 
14.6 ks using ACIS-I and ACIS-S2, with the ACIS-I field 
of view centered at a = 19h21ml2.00s, 6 = -h43°56'43.7", 
roughly on the surface brightness peak. The pixel scale 
is 0".492. Time bins were checked for periods with count 
rates greater or less than 20% of the mean; no such inter- 
vals were found. Hence all of the data with grades of 0, 2, 
3, 4, and 6 were used. The ACIS-I data were adjusted for 
charge-transfer inefficiency (CTI) using the PSU CTI cor- 
rector (Townsley et al. 2000). We used the Chandra data 
analysis software CIAO, version 2.2, for data reduction. 
All spectral analysis was done using the X-ray spectral 
fitting package XSPEC, version 11.2. 
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Fig. 2. — Raw counts ACIS-I image in the 0.5-5 keV band, with 
pixels binned by 4 (i.e., the pixel scale is ~ 2"). The merger feature 
is visible to the southeast of the brightness peak, and the "tail" of 
diffuse emission is seen extending to the northwest. 

2.1. Background and Imaging 

Because A2319 is a large, nearby cluster, its emission 
fills the ACIS-I chip, preventing a direct background mea- 
surement from that data set. The count rate in the S2 chip 
is found to be roughly two times higher than the typical 
background rate, making its use for background estimation 
likewise dubious. One source of this higher than expected 
rate could be a flare affecting the entire observation; how- 
ever, the uniformity of the count rate over the exposure 
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time makes this unlikely, and a visual inspection of the S2 
spectrum does not reveal any flare- like features. A clear 
brightness gradient is visible in the exposure-weighted S2 
image, as well as in the ROSAT PSPC image shown with 
the Chandra footprint in Figure 1; hence, it is clear that 
emission from the very extended cluster is present in the 
S2 data. 

Because there were no portions of the observation with- 
out significant cluster contamination, we use the Marke- 
vitc;li blank-sky observations^. The background was scaled 
up by ~10% after visual comparison of the S2 spectrum 
and the blank sky spectrum, under the assumption that 
emission in the 7-10 keV band is background dominated. 
The recommended procedure for using these blank-sky files 
is to compare the emission in the 10 12 keV band; how- 
ever, the spectral shapes of the S2 spectrum and the back- 
ground spectrum are somewhat better matched in the 7.0- 
10 keV band, and matching the two spectra in the higher 
band results in obvious oversubtraction at energies below 
10 keV. We compared the blank-sky corrected mean sur- 
face brightness in the S2 data to that of the background 
corrected PSPC observation; the Chandra measured sur- 
face brightness is brighter by a factor of ^ 1.5. 

The raw ACIS-I exposure-weighted counts image in the 
0.5-5.0 keV band is shown in Figure 2. The presumed 
merger feature is visible as a sharp surface brightness 
discontinuity to the southeast of the brightness peak. 
The presence of the nicirgcir signatiirc is much clearer 
than in previous X-ray observations; the arclike discon- 
tinuity and the "tail" of emission towards the northwest 
strongly resemble similar features in merging clusters such 
as A2142 and A3667. This is not the possible merger in 
the northeast-southwest direction discussed by, e.g., Fer- 
etti et al. (1997), as it clearly indicates gas movements 
along the axis connecting A2319A and A2319B. 

2.2. Spectral Analysis 

All spectra are fitted using a single-temperature 
MEKAL model, plus components for absorption along the 
line-of-sight and for absorption due to molecular contami- 
nation of the ACIS detector. We fit spectra in the energy 
range 0.9-10.0 keV; poor understanding of the low-energy 
response of ACIS-I prevents us from using data at lower 
energies. 

We first fit for temperature and abundance, fixing the 
hydrogen column density at the Dickey & Lockman (1990) 
value of 8.33 x 10^" cm^^. Fitting over the entire cluster, 
excluding point sources, gives Tx = 11.8±0.2 keV and 
Z = 0.19 ±0.03 (all abundances are in units of solar abun- 
dance; all fitted uncertainties are at the 1 a level), with 
= 1017 for 594 degrees of freedom. This temperature 
is several standard deviations above previously published 
estimates, e.g., Tx = 9.2 ± 0.7 keV determined by Marke- 
vitch et al. (1998) using ASCA data. This spectrum is 
plotted with residuals in Figure 3. 

Previous studies of A2319 have used hydrogen column 
densities in the range (7.85 — 8.9) x 10^° cm~^; often 
the value of A^h used is not provided. By fitting the 
entire cluster spectrum with varying values of A^h, we 
have found that the emission-weighted mean temperature 
varies roughly linearly with A^Hj with the temperature de- 
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creasing by approximately 0.5 keV per 10^" cm~^ (cluster 
temperature uncertainties are generally ~0.2 keV). Fit- 
ting for the column density along with the other param- 
eters yields Tx = 10.6 ± 0.3 keV, Z = 0.20 ± 0.03, and 
A^H = (10.7 ± 0.5) X 10^0 cm-2, with = 999 for 593 
degrees of freedom. 

The Dickey & Lockman (1990) value for the hydrogen 
column density of 8.33 x 10^° cm~^, as well as other val- 
ues used in previous studies of A2319, fall a few standard 
deviations below the range of our fit value. However, un- 
certainties are not readily available for the Hi survey data 
of Dickey & Lockman (1990); moreover, measured A^h val- 
ues in the region of the sky around A2319 vary to levels 
above our fit value. A2319 lies at a fairly low galactic lati- 
tude where there is a significant amount of ISM along the 
line of sight, and the optically thin assumption for deriv- 
ing A^H likely underestimates the true column density by 
a factor of 1.1 1.3 (Dickey & Lockman 1990). Further, 
with a column density this high there is likely to be a 
significant contribution (> 10%) to the hydrogen column 
by molecular hydrogen (Lockman 2003). Also, fitting A^h 
along with other parameters in our temperature mapping 
suggests that there may be a gradient with magnitude of 
a few 10^" cm~^ across the ACIS-I image. For the rest 
of the paper we adopt the value of 8.33 x 10^° cm~^, but 
readers should keep in mind that it is almost certainly too 
low. 
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Fig. 3. — Entire cluster spectrum (excluding point sources) and 
residuals, plotted with the best-fit MEKAL spectrum described in 

the text. 

The uncertainty of 0.2 keV given for the cluster temper- 
ature above includes only the statistical uncertainty from 

the spectral fit. We adopt a l-cr uncertainty in A'h of 
~ 10^° cm~^, which introduces a corresponding 0.5 keV 
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uncertainty in the temperature. The background subtrac- 
tion also affects the temperature. The Poisson uncertainty 
in the background scale factor determined using the 7~ 
10 keV S2 spectrum is ^ 4%, corresponding to a 0.2 keV 
uncertainty in the cluster temperature. In addition, the 
background scaling using the 7-10 keV band produces a 
cluster temperature that is 0.3 keV higher than that when 
scaling the background using the 10-12 keV band; thus, 
we adopt a temperature uncertainty contribution from 
the background scaling of 0.3 keV. Combining our three 
sources of uncertainty (statistical, Nn, and background 
scaling), we arrive at a cluster temperature and uncer- 
tainty of 11.8±0.6 keV. It should be noted that hydrogen 
column density uncertainties are not included in tempera- 
ture uncertainties in the rest of the paper unless explicitly 
noted. 

Because X-ray point sources are visible in the Chandra 
data that were not noticeable in previous observations, it 
is possible that their presence could have affected previous 
temperature measurements. To check this, we also fit the 
entire cluster spectrum without removing point sources; 
this produces a temperature decrease of less than 0.1 keV. 

Our measured value of Tx = 11.8 ± 0.6 keV is some- 
what higher than previous temperature measurements; our 
abundance value of 0.19±0.03 is low in comparison to pre- 
vious studies, though abundances in this range appear to 
be typical in studies of merging clusters (De Grandi & 
Molendi 2001). The discrepancy between our temperature 
measurement and previously published temperatures may 
be partially explained by Chandra^s relatively small field of 
view and the large angular extension of A2319. As is clear 
from Figure 1, there is significant cluster emission outside 
of the ACIS-I field. Using the PSPC image, we found that 
^ 30% of the cluster emission in the 0.5-2.0 keV band 
lies outside of our ACIS-I observation. A MEKAL model 
fit on the S2 chip (excluding point sources) gives a tem- 
perature of 7.1 ± 1.2 keV (x^ = 214 for 204 d.o.f.). This 
value is in agreement with ASCA measurements of 6-9 




keV in large regions around and including the area cov- 
ered by our S2 observation (Markevitch 1996). If the bulk 
of the gas outside the ACIS-I field is similarly cooler than 
our measured average temperature of the cluster, then our 
temperature measurement with Chandra would naturally 
be higher than measurements with previous-generation, 
larger field of view instruments. This effect probably does 
not account for the entire difference between our result 
and others, because measurement of temperatures within 
small regions of the ACIS-I chip give slightly higher-than- 
expected results as well, as will be discussed in §2.3. 

If a higher value for TVh were used, as discussed above, 
our fit temperature would be lower. This cannot ac- 
count for the discrepancy between our results and previ- 
ously published measurements, however, as previous stud- 
ies have used column densities within ^ 0.5 x 10^" cm~^ 
of our adopted value. 

The cluster temperature fit is sensitive to the choice 
of energy band. For example, fitting the entire cluster 
spectrum (with abundance and hydrogen column den- 
sity allowed to vary) in the 0.9-10.0 keV band gives 
Tx = 10.6 ± 0.3 keV (x^ = 999 for 593 d.o.f.); however, 
fitting between 1.7-10.0 keV gives Tx = 6.2 ± 0.2 keV 
(;^2 ^ 705 for 559 d.o.f.), and fitting between 2.0-10.0 keV 
gives Tx = 7.6 ± 0.4 keV (x^ = 576 for 523 d.o.f.). While 
the specific behavior will vary by instrument, it should be 
noted that the lower energy limit of most previous temper- 
ature measurements has been ~1. 5-2.0 keV. One obvious 
explanation for the extreme dependence of spectral fitting 
on energy band is simply that the cluster is not isothermal, 
as we show in §2.3. 

2.3. Temperature Structure 

Chandra provides the means to perform a much more 
detailed study of the temperature structure of A2319 than 
previous instruments, permitting inspection of the cluster 
merger features. To this end we have made an X-ray tem- 
perature map of A2319 using the ACIS-I data. The map 




Fig. 4. — X-ray temperature map (left), and significance map (right). The contours are from the 0.5—5 keV energy band image shown in 
Figure 2 after smoothing with a Gaussian of constant size, and are spaced at 10% of peak cluster intensity. Temperature pixels are 1' (66 
kpc) on a side. The average temperature (T) = 11.8 keV, and uncertainties in this average temperature are not included in the significance 
map. 
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was created by measuring the temperature at each point 
on a grid, using a circular region enlarged until it con- 
tained 2000 counts in the 0.9-2.0 keV energy range. The 
regions overlap slightly at the center, and increasingly to- 
wards the edge; hence the pixels are not independent of 
one another. In the faint regions of the observation, where 
fitting region radii are larger than two pixel widths, only 
one pixel in four is measured. The spectra at each point 
were calculated using the same procedure as for the whole 
cluster spectrum described in §2.2, with abundance float- 
ing and Nu = 8.33 x 10^° cm~^. The abundance was left 
as a free parameter as abundances are known to vary in 
merging systems; fixing it to the cluster average produces 
temperature changes of < la across the temperature map. 




Fic;. 5. — V-band image from the Digitized Sky Survey with Chan- 
dra observation overlaid. Contours are the same as in Figure 4. The 
central (i.e., brightest) galaxy of A2319B is indicated with an arrow. 
The bar indicates a distance of 0.5 Mpc. 

The temperature map is shown in Figure 4 (left) with 
overlaid surface brightness contours. Also in Figure 4 
(right) is a map of the significance of deviations from the 
mean temperature; that is, the difference between each 
pixel temperature and our adopted cluster mean tem- 
perature of 11.8 keV, divided by the uncertainty in the 
pixel temperature. The general structure of the tempera- 
ture map includes two coolcr-than-avcragc regions that lie 
along a northwest-southeast line, and possibly two hotter- 
than-average regions that lie to the northeast and south- 
west of center. This temperature morphology is suggestive 
of a merger along a northwest-southeast trajectory, where 
remnants of cold cores remain and shock heated gas is 
escaping perpendicular to this merger axis, as seen in hy- 
drodynamical simulations (Roettiger et al. 1997; Ricker & 
Sarazin 2001). The cold spots deviate from the mean by 
> 2a; the hot spots are somewhat less significant. The 
very high 15 keV) temperature regions lie where the 
cluster surface brightness is lowest, making these temper- 
atures particularly susceptible to background subtraction 
errors. Overall, temperatures are higher than would be 
expected based on previous studies of A2319 (Markevitch 
1996). Regardless of any overall temperature increase, the 
nonisothermality of the cluster provides some indication 
as to the origin of the poor fit discussed in §2.2. 



The level of substructure revealed here is more detailed 
than has been previously seen. The coolest region lies just 
south of the surface brightness peak, perhaps indicating a 
cool core that has thus far survived the ongoing merger. It 
is not immediately obvious from this map whether there is 
a sharp temperature change across the merger feature sig- 
nificant enough to deduce the existence of either a shock 
front or a cold front; we examine this in more detail in 
§3.1. 

This cool core; has not been identified in the earlier 
ASCA temperature map (Markevitch 1996). It seems 
likely that surrounding areas of higher-than-average tem- 
peratures obscured the core in the lower angular resolution 
ASCA map. Molendi et al. (1999) pointed out a "subclus- 
ter" of temperature 6.9 ± 1.0 keV to the northwest of the 
cluster center, and there is evidence for the presence of 
this cool region in the temperature map of Markevitch 
(1996). This subcluster is seen here ^ 6' northwest of the 
X-ray brightness peak, though at a somewhat higher tem- 
perature. Also present is a distinct region of somewhat 
elevated (i.e., above the mean) temperatures between this 
subcluster and the cool center. 

The cool ICM "subcluster" has been identified with sub- 
cluster A2319B; however, at this resolution it is clear that 
the cool region is not associated with the center of A2319B, 
but rather lies 2-3' to the east-southeast of it, as can be 
seen by comparing the temperature map to the visual- 
band image shown in Figure 5. This suggests that the 
subcluster is in a transient phase wherein the ICM has de- 
coupled from the galaxies. Such a state has been observed 
in other merging clusters such as 1E0657-56 (Markevitch 
et al. 2002) A754 (Zabludoff & Zaritsky 1995; Markevitch 
et al. 2003), CI J0152.7 1357 (Maughan et al. 2003), and 
A2034 (Kempner et al. 2003). 

Overall, the temperature map reveals complex substruc- 
ture of the type now known to occur in galaxy cluster 
mergers. Such substructure is also seen in hydrodynam- 
ical simulations (e.g. Roettiger et al. 1997; Onuora et al. 
2003) 

3. MERGER ANALYSIS 

Wc present here a simple analysis of the merger features 
in A2319, wherein we assume a simple spheroidal geom- 
etry for an isothermal body of gas falling into a relaxed 
/3-model cluster. This is what might be termed the "tradi- 
tional" analysis of a merging cluster (following Vikhlinin 
et al. 2001). However, numerical simulations of clusters 
(e.g., Ricker & Sarazin 2001; Bialek et al. 2002; Nagai 
& Kravtsov 2003; Onuora ct al. 2003) have made it clear 
that the dynamics within a mid-stage merger are much 
more complex than this. Nevertheless, this naive analysis 
offers some level of quantitative information about the na- 
ture of the merger front, and permits comparison to other 
merger analyses. 

3.1. Temperature and Brightness Profiles Across Merger 

Feature 

We measure the surface brightness and temperature pro- 
files across the merger feature (see Figure 6). The bright- 
ness is measured in arcs on a wedge, chosen with a radius 
of curvature and angular width that match the brightness 
discontinuity reasonably well. We then measure the tern- 
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perature, making spectra as previously described, in arc 
segments of sizes chosen both to provide a sufficient num- 
ber of photons and to permit study of temperature vari- 
ation across the front; we select the segment boundaries 
to avoid having a region straddling the surface brightness 
discontinuity. Note that this is not a cluster radial profile; 
the wedge in which this is measured is chosen to match 
the brightness discontinuity, and is not centered on the 
brightness peak. 

While there is clearly a brightness change, this change 
is not as sharp as those seen in merging clusters such as 
A3667 (Vikhlinin et al. 2001). This can be readily ex- 
plained if the merger is not taking place close to perpendic- 
ular to the line of sight; indeed, the aforementioned differ- 
ence in line-of-sight velocity between A2319A and A2319B 
of ~ 2900 km s~^ (Oegerle et al. 1995) suggests that we are 
viewing the merger at some large angle. This introduces 
substantial uncertainties into the analysis below. 
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Fig. 6. — Brightness (in arbitrary units) and temperature profiles 
across the merger feature. The arrow indicates the approximate po- 
sition of the brightness discontinuity seen in Figure 2. A fit to the 
surface brightness inside the merger feature using Eq.2 is shown by 
the sohd hne. 

If we assume that the infalling gas body is a spheroid 
with constant gas density, then the surface brightness pro- 
file at distances from the front much less than the major 
axis of the spheroid will be given by 



S{d) = 2^/'^\/Reoy/d , 



(1) 



where R is the radius of curvature at the front, eq is the 
volume emissivity of the gas, and d is the distance from the 
front (Vikhlinin et al. 2001). This function adequately de- 
scribes the surface brightness profile of our data in the re- 
gion just inside the front (i.e., the region between roughly 
1.5-2.8' in Figure 6). 
More precisely, the surface brightness profile is 
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for \f3\ < 0.25, where a and b arc the short and long axes 
of the spheroid, respectively (Vikhlinin et al. 2001). Fit- 
ting this function to the brightness profile just interior to 
the front gives /3 < 0.1. Our approximation of constant 
density inside the front is thus justified. 

An examination of Figure 6 docs not conclusively deter- 
mine the nature of the merger feature, i.e., whether it is 
a shock front or cold front. The temperature just inside 
the merger is 9.0 ± 0.9 keV, while the temperature just 
outside is 10.8 ± 1.3 keV. This suggests that the feature is 
a cold front, but the temperature uncertainties are large. 
The temperature falls by another 1-2 keV deeper inside 
the infalling subcluster, perhaps suggesting some heating 
near the surface brightness discontinuity. 

To rule out significant upward biasing of temperatures 
inside the front by projected hotter gas in front of and 
behind the cooler gas, we fit a two-component MEKAL 
model to a region inside the front, near the brightness dis- 
continuity, with the hotter component fixed to the temper- 
ature measured outside the discontinuity. We find that to 
measure a cool component temperature that is 1 cr lower 
than the single-component temperature measurement re- 
quires a hot component contribution of > 40% of the emis- 
sion. As this seems unreasonably high, we conclude that 
our temperature measurements inside the front are not 
significantly biased by projected hotter gas. 

3.2. Density Variation Across Merger Feature 

In general, the intensity of a body of gas at constant 
temperature is 



1 



■ J nenHA{Tx,l)dl , (3) 



where A{Tx,l) is the emissivity of the gas and the length 
element dl is along the line of sight; the integration is car- 
ried out over the entire body along the line of sight. 

If the spheroid's long axis is much larger than the minor 
axes, we can model the infalling subcluster as a "bullet" of 
width L; we assume a constant temperature and intensity. 
Using these assumptions in equation (3) and solving for 
the electron density gives 



Me 



LA(rx) 



1/2 



(1+-)^ 



(4) 



We use fie = 1-67 and fiH = 1-4, the values for a fully 
ionized gas of one-third solar abundance. Using estimated 
values for / and L, we obtain an electron density immedi- 
ately inside the merger front of (6.0 ± 1.0) x 10~^ cm~^. 

To get the electron density outside the front, we assume 
that the gas fits a spherical /3-model density profile: 




-3/3/2 



(5) 



with central electron density rieo ^^'^ critical radius 9c- 
That is, we assume that the gas to the southeast of the 
front is part of the original relaxed cluster into which a 
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subcluster is falling, and is thus far unperturbed by the 
merger. To get values for (3 and 9c, we fit the surface 
brightness in our wedge, outside the merger. We arrive at 
an electron density immediately outside the merger fea- 
ture of (2.0 ± 0.6) X 10^'^ cm^"', or roughly one-third the 
density immediately inside the feature. 

These densities correspond to pressures {p = ngTx) in- 
side of Pin = (5.4 ± 1.0) X 10~^ keV cm~^ and outside of 
p^nt = (2.2±0.7)x 10^^ kcV cm^-^. Using the relationships 
between these pressures and the Mach number M of the 
infalling gas cloud (where M = v/cout is the Mach number 
in the free stream outside the merger) gives M = 1.1 ± 0.3 
(sec §122; Landau & Lifshitz 1987). The infalling cluster 
would thus probably be moving at a roughly sonic speed 
if indeed the merger were taking place in the plane of the 
sky, as we have assumed for this analysis. Because of the 
line of sight velocity difference of the galaxies associated 
with A2319A and A2319B, wc expect that the merger axis 
does not lie in the plane of the sky. 

3.3. Toward a Cluster Dynamical Model 

Combining the results of the previous sections, we 
present the following picture of the merger in A2319. 
There is a jump by a factor of 3.0 ± 1.0 in the density 
of the gas as one crosses the brightness discontinuity from 
the unperturbed gas outside the merger towards the clus- 
ter core. This is accompanied by a slight temperature 
decrease of 1-3 kcV, and a brightness increase by a factor 
of ^ 3; the combined densities and temperatures give a 
pressure jump by a factor of 2.5 ± 0.9. These results indi- 
cate the presence of a cold front, although the temperature 
difference across the front is not as large as is observed in, 
e.g., A3667 (Vikhlinin et al. 2001). 

However, we have assumed for this analysis that the in- 
falling subcluster is moving in the plane of the sky; our 
value for the electron density in the cool core is thus an 
overestimate given the known line-of-sight velocity differ- 
ence between A2319A and A2319B that indicates that 
bulk gas motions are not perpendicular to the line of 
sight. This is most easily seen by examining equation 4; 
if the subcluster is not oriented perpendicular to the line 
of sight, then we are overestimating the X-ray intensity I, 
and hence also the electron density rig. Moreover, if this 
merger has a nonzero impact parameter, then our estimate 
for the ambient electron density, i.e., the density outside 
the merger feature, is likewise an overestimate. It is thus 
possible that the inside/outside density and pressure ra- 
tios are in fact lower or, less likely, higher than the values 
given. This does not change the general interpretation of 
the merger feature as a cold front; the temperature change 
is indisputable, and the uncertainties in density are not 
large enough to accomodate a pressure outside the front 
greater than that inside the front. 

The simplest interpretation for the merger geometry 
seen in A2319 is that A2319B has recently fallen through 
A2319A, in the process losing much of its ICM as indi- 
cated by the low X-ray brightness around the giant ellipti- 
cal that dominates its galaxy population. The separation 
of the cold spot near A2319B from its galaxies supports 
this (see Figures 4 & 5). However, the structure of the cold 
front suggests motion away from A2319B. We suggest that 
the encounter of the two subclusters of A2319 has caused 



the cool core of A2319A to be displaced from its position 
at the center of the subcluster, and that this core is now 
recoiling from that displacement and has passed its orig- 
inal, central position. This is supported by the fact that 
the coldest part of A2319A's core is located slightly to the 
southeast of the brightest cluster galaxy. The merger fea- 
ture is then a result of the interaction of the dense core 
ICM with less dense, warmer ICM surrounding the core. 

The apparent survival in some form of the cold core 
ICM of A2319B may indicate a non-zero impact parame- 
ter. Given this and the relative sizes of the two subclusters, 
displacement of the core of A2319A to the point of creat- 
ing motion of the core at near-sonic speeds would require 
a large infall velocity. It is also possible that the merger 
was essentially head-on, and that cool ICM spatially asso- 
ciated with the galaxies of A2319B is not from the original 
core of the subcluster, but was pulled from the core ICM 
of A2319A during the collision (see Pearce et al. 1994). 

We estimate the timescale since closest approach of the 
two subclusters by constructing a simple, two-body dy- 
namical model. Using the line of sight velocity dispersion 
of A2319A {aA = 1324 km s-\ Oegerle et al. 1995), we 
obtain a crude estimate of the collision infall velocity of 
•\/6c7A = 3243 km s~^ (this assumes infall from infinity). 
Combined with the measured line-of-sight velocity differ- 
ence of subclusters A and B (2909 km s~^ Oegerle et al. 
1995), we estimate that the merger trajectory has an angle 
of ^ 65° out of the plane of the sky. The corresponding 
velocity in the plane of the sky is ~ 1430 km s~^. This 
gives a time since closest approach of the subcluster cores 
of - 0.4 Gyr. 

We emphasize that this is only one possible merger sce- 
nario. It does not include the possibility of a second merger 
event taking place along the northeast-southwest direction 
such as that suggested by an analysis of earlier X-ray data 
(Feretti et al. 1997). 

4. CLUSTER OBSERVABLES DURING A MAJOR MERGER 

The merger signatures in Abell 2319 are clear. These 

include significant centroid shifting in the Einstein IPC 
(Mohr et al. 1995) and ROSAT PSPC X-ray images; two 
subclusters identified in the optical (Faber & Dressier 
1977; Oegerle et al. 1995); differing distributions of galax- 
ies and ICM; and a surface brightness discontinuity and 
temperature structure in the Chandra data. We seek now 
to examine how mergers perturb the global physical struc- 
ture of clusters. We address this empirically by simply 
examining how particular bulk properties of A2319 (bind- 
ing mass, ICM mass, isophotal size, luminosity, emission 
weighted mean temperature, and galaxy light) compare 
to typical galaxy clusters. Specifically, we compare the 
properties of A2319 to what is essentially a flux-limited 
ensemble of 44 galaxy clusters from the nearby universe 
(Mohr et al. 1999). While the high resolution of Chandra 
is not necessary for this, the Chandra observation nonethe- 
less provides another high-quality data set for such study. 

The question of how much merging perturbs the global 
structure of galaxy clusters is particularly important in 
light of the planned and ongoing high-yield galaxy cluster 
surveys. In these surveys, rather simple observables like 
the SZE flux. X-ray flux, and galaxy light will be used to 
estimate cluster masses for studies of the dark energy (e.g., 
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Haiman et al. 2001). Even though it has recently been 
shown that very large surveys contain enough information 
to self-calibrate while precisely constraining the dark en- 
ergy (Majumdar & Mohr 2003a, b; Hu 2003), any improve- 
ments in our understanding of cluster mass observable re- 
lations, their evolution, and the effects of merging on them 
will lead to tighter limits on systematic uncertainties in the 
resulting cosmological constraints. 



radius [Mpc] 
0.4 



> 




radius [arcmin] 

Fig. 7. — Projected temperature profile created around the cen- 
ter of the cluster as found by a /3-modcl fit, not around the surface 
brightness peak. The outermost annulus partially extends beyond 
the boundaries of the ACIS-I observation, and so we mark the an- 
gular extent of this annulus with a dashed line. 



4.1. Naive Analysis of Mass Profile 

We obtain a naive measure of M2500, i-e., the mass en- 
closed by r25oo, the radius within which the mean density 
is 2500 times the critical density of the universe. Most clus- 
ter surveys focus on properties at larger radii (e.g., rsoo), 
but Chandra^s small field of view relative to the large an- 
gular extent of A2319 makes this impossible with a single 
pointing. 

We assume that the cluster ICM density distribution is 
fit by a spherical /3-model and that the ICM is in hydro- 
static equilibrium; the binding mass within a radius r is 
then 



jinipG \ dhir 91nr / 

Figure 7 contains a projected temperature profile of 
A2319, where the cluster center is that found by the (5- 
model fit to the X-ray surface brightness (described be- 
low). There is no easily quantifiable variation of tem- 
perature as a function of distance from the cluster cen- 
ter. Given the measured density and temperature pro- 
file, r"25oo lies mostly outside the ACIS-I image. To es- 
timate the mass at this radius, we adopt an isothermal 
temperature profile and extract the average temperature 
from the outer three annuli in Figure 7. This temperature 
is 11.1 ± 0.9 keV, less than la lower than the emission 
weighted mean temperature for the cluster. 



A fit to the Chandra surface brightness image of A2319 
gives core radius r^ = 0.17 ± 0.01 Mpc {9^ = 2'.6 ± O'.l) 
and /3 = 0.55±0.01 (compared to fc = 0.15±0.05 Mpc and 
f3 = 0.54±0.06 from the analysis of the PSPC image; Mohr 
ct al. 1999). With these values and equation (6), we find 
r2500 = 0.67 ±0.02 Mpc (6*2500 = 10'.2±0'.4) and binding 
mass M2500 = (4.2 ± 0.5) x 10^^ Mq. The uncertainties 
quoted here for the /3-model fit are 1 a statistical uncer- 
tainties only, and do not reflect the fact that the f3 model 
is not a particularly good fit the surface brightness in this 
complex cluster. The mass uncertainty is dominated by 
the uncertainty in the temperature measurement. 

We compare this mass estimate to that expected for a 
cluster with this emission-weighted mean temperature, us- 
ing an Af25oo-2x relation derived from a sample of seven 
intermediate-redshift clusters (Allen et al. 2001). For our 
cluster temperature of Tx = 11.8 ± 0.6 keV the best fit 
relation gives M2500 = (6.7 ± 0.8) x lO^'^M©, which is a 
factor of 1.6 ± 0.3 higher than our value. In this merging 
cluster, the hydrostatic equilibrium assumption and spher- 
ical /3-model fitting thus lead to a mass estimate that lies 
~ 60% off the relation found in apparently "relaxed" clus- 
ters. The Allen ct al. (2001) sample is too small to make 
meaningful statements about the scatter, but other anal- 
yses of much larger samples show scatter at roughly the 
25% level (Finoguenov et al. 2001). 

4.2. Comparison of A 23 19 to Large Cluster Sample 

We examine five bulk properties of A2319 — the X-ray 
luminosity, emission-weighted mean temperature, ICM 
mass, isophotal size, and K-h&nd galaxy light — and com- 
pare these properties to the same properties for large sam- 
ples of galaxy clusters. In the case of all but the galaxy 
light we use an ensemble of 44 clusters studied using the 
ROSAT PSPC (Mohr & Evrard 1997; Mohr et al. 1999, 
2000), but reanalyzed at the cluster radius r25oo- For each 
cluster wc determine r25oo using the emission weighted 
mean temperature and the published M25oo-temperature 
relation (Allen et al. 2001). For the PSPC sample, expo- 
sure corrected, background subtracted images were pre- 
pared in the rest frame 0.5-2.0 keV band for each cluster. 
In the case of the galaxy light, we compare to an ensemble 
of ~100 clusters whose properties are being studied using 
X-ray data and 2MASS near-IR data (Lin et al. 2003). 

We measure ^2500, the luminosity projected within a cir- 
cle of radius ^2500 in the 0.5-2.0 keV band. Due to the un- 
certainty of the spectral response of Chandra below ~ 0.9 
keV, for our observation of A2319 we measure the flux or 
luminosity in an image that includes only counts in the 
0.9-2.0 keV band. With the emission-weighted mean tem- 
perature, we calculate the conversion between the count 
rate in this band and the flux within the rest frame 0.5-2.0 
keV band. Another difficulty is that we find 62500 = ll'-8 
from the Allen et al. (2001) A'hcMO-Tx relation, slightly too 
large to fit within the ACIS-I observation. However, the 
low luminosity near the edges, relative to the central lu- 
minosity, means that our value Z/2500 = 5.2 x 10*^ erg s~^ 
contains the bulk of iv2500- Indeed, using the Chandra foot- 
print on the i?OS'^TPSPC image, we find that 15% of the 
flux within r25oo is missed; thus, our corrected estimate of 
the luminosity in the 0.5-2.0 keV band is 1/2500 = 6.0 x 10^'' 
erg s~^. This is high by ^18% compared to the value 
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^2500 = 5.1 X lO'*'^ erg measured using the PSPC im- 
age. Figure 8 (top) contains the ROSAT sample (small 
points) with best fit power law together with the Chandra 
measurement (large point). The luminosity is low by 54% 
relative to the expectation for a cluster with a temperature 
of 11.8 keV, compared to an RMS fractional scatter about 
the best fit relation of 57% (the PSPC value for L2500 is 
low by - 61%). 
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Fig. 8.— Scaling relations for the sample of ROSAT PSPC ob- 
servations from Mohr et al. (1999) (small points), plus our mea- 
surements for A2319 {Chandra large, solid point; PSPC large, open 
point), using a temperature of 11.8 keV. Best fits to the PSPC sam- 
ple are shown as lines. Top: L2500 in the 0.5—2.0 keV band. Mid- 
dle: MjcM within r2500. Bottom: Isophotal radius for an isophote 
of 1.53 X 10"^'^ erg s~^ cm~^ arcmin"-^ in the 0.5-2.0 keV band. 
Points shown in outline were excluded from the fit, as the use of a 
very high isophote caused them to give erroneous results. 



We measure the ICM mass within r25oo using a mea- 
surement of the flux from the cluster combined with the 
/3-model fit parameters and our cluster temperature of 
11.8 keV. The count rate emissivity of a parcel of gas 



within the 0.5-2.0 keV band has low sensitivity to tem- 
perature variations and, assuming all the ICM is emitting 
at the emission weighted mean temperature, provides a 
good estimate of the ICM mass (see Fabricant et al. 1980; 
Mohr et al. 1999). The ICM mass from the Chandra anal- 
ysis is MicM = 6.9 X IO^^Mq, corresponding to an ICM 
mass fraction of /icm — 16%. The corresponding value 
from the PSPC analysis is AficM = 5.2 x lO^^Af©, roughly 
30% lower. Figure 8 (middle) contains the i?05'j4T' sample 
(small points) with best fit power law together with the 
Chandra measurement (large point). The ICM mass is low 
by ~ 14% relative to the expectation for a cluster with a 
temperature of 11.8 keV, compared to an RMS fractional 
scatter about the best fit relation of 20% (the PSPC value 
for the ICM mass is low by ~ 35%). 

We measure the isophotal size Ri for A2319 at an 
isophote of 1.53 x 10"^'^ erg s~^ cm~^ arcmin^^ in the 
0.5-2.0 keV band. This isophote is chosen so that the 
isophotal size is not affected by the limited field of view of 
the Chandra footprint. We measure the size using the area 
Al enclosed by this isophote, and find an equivalent radius 
from Al = 7ri?j. For the Chandra observation we obtain 
Ri = 0.44 Mpc, compared to the value Ri = 0.39 Mpc 
obtained when using the PSPC image. Figure 8 (bottom) 
contains the ROSAT sample (small points) with best-fit 
power law together with the Chandra measurement (large 
point). Note that at an isophote this bright there are sev- 
eral clusters that simply fall off the relation defined by the 
bulk of the PSPC clusters. We are forced to use such a 
bright isophote because the Chandra footprint is so small 
compared to the angular extent of A2319. Nevertheless, 
A2319's isophotal size is ~ 14% lower than that expected 
for a cluster with its emission weighted mean temperature, 
compared to an RMS fractional scatter about the best fit 
relation of 14% (the PSPC value for Ri is low by ~ 24%). 

Interestingly, A2319 does not stand out significantly 
from the sample of 44 clusters (essentially an X-ray fiux- 
fimited sample) studied with the ROSAT PSPC. The 
merger in A2319 does not perturb the cluster significantly 
in luminosity, ICM mass, or isophotal size from the val- 
ues expected for a cluster with its emission weighted mean 
temperature. In addition, an analysis of the galaxy light in 
the if-band that is projected within rsoo in A2319 leads to 
an estimate of the cluster X-band light that is 14% higher 
than expected for a cluster with a 11.8 keV temperature, 
when compared to a sample of ^ 100 clusters where the 
rnis scatter is 30% (Lin et al. 2003). 

One possible explanation is that the merger event is rel- 
atively minor (the ratio of velocity dispersions of A2319A 
and A2319B suggests a mass ratio of ^8), but it may 
also be that merging clusters are perturbed in all their 
quantities in such a way that they remain close to the 
population-defined scaling relations. In fact, it should be 
noted that many of the clusters contained in the PSPC 
sample exhibit evidence for ongoing mergers (Mohr et al. 
1995; Buote & Tsai 1996). We cannot hope to deliver a fi- 
nal verdict on the effects of merging on the bulk properties 
of clusters with studies of single clusters; however, our re- 
sults do provide some evidence that bulk properties either 
do not change much as a result of mergers, or change in a 
correlated way that maintains the strikingly small scatter 
of scaling relations. 
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Correlated changes in luminosity and temperature 
within merging clusters have been examined with nu- 
merical simulations. Ricker & Sarazin (2001) measured 
the luminosity and temperature boosts in merging clus- 
ter systems as a function of time. If we assume a 1:3 
mass ratio for the subclusters, then the simulations pre- 
dict a peak luminosity boost by a factor of ~ 2-4, along 
with a peak temperature boost of a factor of ~ 1.5-2.0, 
with correspondingly smaller boosts associated with larger 
mass ratio mcrgcirs. Simultaneous boosts to the lumi- 
nosity and temperature of these magnitudes would not 
make A2319 stand out in the luminosity-temperature re- 
lation in Figure 8. However, A2319 appears "normal" 
with respect to its luminosity-temperature, isophotal size- 
temperature, ICM mass temperature and galaxy light- 
temperature properties. It would seem to be contrived 
to claim that large, merger related excursions in these five 
bulk properties of the cluster all take place in just such 
a way as to keep the cluster near the observed, typical 
behavior for a large sample of clusters. A simpler expla- 
nation would appear to be that these five cluster properties 
are simply not dramatically affected by the merger taking 
place in A2319. 

5. CONCLUSIONS 

Using Chandra data, we have identified and studied a 
major merger event in A2319 that appears to be taking 
place along the axis connecting its two major optical sub- 
clusters. The X-ray brightness map shows a clear discon- 
tinuity that appears similar to cold fronts found in other 
clusters. Although this cold front appears to be as large as 
the one studied in A3667 (Vikhlinin et al. 2002), it is not 
as sharp. This, together with previous measurements of 
the line of sight velocity difference between the two main 
optical subclusters (Oegerle et al. 1995), suggests that the 
merger is not taking place in the plane of the sky. We 
propose a merger model where the trajectory lies approxi- 
mately 65° out of the plane of the sky, and at this viewing 
angle it becomes even more challenging to make quanti- 
tative statements about the ICM properties near the cold 
front. Nevertheless, we estimate that the prcssiirc change 
across the front is < 2.5, and that the higher density ICM 
also has the lower temperature. The estimated merger 
Mach number of ^ 1.1 is likewise consistent with other 
merging systems such as A3667. We propose a two body 
merger where A2319B merged from the southeast travel- 
ing northwest, with the A2319B galaxies and dark matter 
passing through the A2319A core roughly 0.4 Gyr ago. 

The measured emission-weighted mean temperature of 
this messy, merging cluster is Tx = 11.8 ±0.6 keV and the 
mean abundance is Z = 0.19 ± 0.03, using a hydrogen col- 
umn toward the cluster of A^h = 8.33 x 10^° cm'^. The fit 
values deviate somewhat from previous studies of A2319. 
Our higher temperature is likely due in part to the; small 
field of view of Chandra compared with other instruments 
used to study this cluster. Also, we have shown that the 



emission-weighted mean temperature depends sensitively 
on the choice of energy band, which is at least partly ex- 
plained by the highly nonuniform temperature structure 
we have revealed in this cluster. 

Our temperature map shows substructure now consid- 
ered typical in merging systems. The cool core of A2319A 
is readily visible, and the angular separation of the galax- 
ies of A2319B and an associated cool ICM region indicates 
a separation of the galaxies from the ICM of this subclus- 
tcr, a transient phenomenon that gives further evidence of 
a merger event. There is some evidence for a hot bridge 
of ICM between the two cores, a characteristic associated 
with shock heating in mergers that has been seen in sim- 
ulations. 

We examine how this merger affects the bulk proper- 
ties of A2319. We naively apply the hydrostatic equi- 
librium assumption to measure a total mass within r25oo 
of M2500 = (4.2 ± 0.5) X 10" Afe, a factor of 1.6 ± 0.3 
lower than the mass predicted by a mass-temperature rela- 
tion derived from five intermediate-redshift clusters (Allen 
et al. 2001). This offset is the strongest indication that 
the structure of Abell 2319 has been significantly affected 
by the merger. Our measurements for L2500 temperature, 
AficM-temperature, and isophotal size-temperature are 
compared to a sample of 44 clusters observed with ROSAT 
PSPC (Mohr & Evrard 1997; Mohr et al. 1999, 2000). In 
all three cases our measured values for A2319 are within 
the scatter of the PSPC-dcrived scaling relations. In ad- 
dition, we note that the if-band light in the A2319 galaxy 
population is consistent with that expected for a cluster of 
this emission weighted mean temperature (Lin et al. 2003). 
It is possible that changes in bulk parameters due to merg- 
ers are actually quite large but take place in a correlated 
way that maintain the low, observed scatter in cluster scal- 
ing relations; this has been shown for some properties in 
numerical simulations (Ricker & Sarazin 2001; Evrard & 
Gioia 2002). However, it will require further studies to de- 
termine whether it is possible for mergers to create large, 
correlated displacements in five cluster parameters (i.e., 
emission- weighted mean temperature, luminosity, isopho- 
tal size, ICM mass, and K-band galaxy light) that main- 
tain the low scatter in all four scaling relations. Another 
possibility is that despite the X-ray imaging spectroscopy 
and optical evidence for an ongoing merger in A2319, a 
merger of this scale is simply not sufficient to grossly per- 
turb the bulk properties of the cluster. 
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Digitized Sky Surveys were produced at the Space Tele- 
scope Science Institute under U.S. Government grant NAG 
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